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ABSTRACT 

We employ a 3-dimensional magnetohydrostatic model of a horizontal flux tube, embedded in a mag- 
netic surrounding atmosphere, to successfully reproduce the azimuthal and center-to-limb variations of 
the Net Circular Polarization observed in sunspot penumbrae. This success is partly due to the realistic 
modeling of the interaction between the flux tube and the surrounding magnetic field. 

Subject headings: Sun: sunspots - Sun: magnetic fields - Sun: polarimetry 



1. INTRODUCTION 

A possible scenario that explains the magnetic field con- 
figuration of the sunspot penumbra is the so-called un- 
combed penumbral model (Solanki & Montavon 1993). In 
this model a horizontal flux tube that harbors the Ever- 
shed flow is embedded in a more vertical magnetic field. 
A similar configuration has been obtained by Heinemann 
et al. (2007), who carried out 3D MHD simulations of the 
penumbra and found that penumbral filaments are pro- 
duced by bubbles of weak and horinzontal magnetic field 
fully embedded in a stronger and more vertical one. They 
also find field free gaps connected to the deeper convec- 
tion zone which appear as bright regions in the emergent 
intensity. 

One of the most critical observations that any penum- 
bral model must reproduce, is the Net Circular Polar- 
ization (Sanchez Almeida & Lites 1992). Different re- 
alizations of the uncombed penumbra have successfully 
explained these observations (Solanki & Montavon 1993; 
Martinez Pillet 2000; Schlichenmaier et al. 2002; Miiller et 
al. 2002, 2006). However, Spruit & Scharmer (2006) have 
pointed out that the models of the uncombed penumbra 
published so far, do not consider the effects of the exter- 
nal field wrapping around the flux tube. The goal of this 
work is to remove this point of inconsistency. As it will be 
shown, a uncombed model that consistently considers the 
perturbation in the external field introduced by a cylin- 
drical flux tube, is also able to explain the Net Circular 
Polarization observed in the sunspot penumbra with sur- 
prising accuracy. 

2. FLUX TUBE MODEL AND REFERENCE FRAME 

Our embedded flux tube model is adopted from Borrero 
(2007). This model considers a cylindrical flux tube that 
carries the Evershed flow. The flux tube is located per- 
pendicularly to the vertical z-axis and is embedded in an 



inclined potential magnetic field. 

This model prescribes the magnetic field configuration 
in the local reference frame (LRF): S = {ex,ey,ez}; see 
Fig. 1), being the flux tube's axis oriented along e^. Ana- 
lytical expressions for Bx{y,z), By{y^z), and Bz{y^z) are 
taken from Eqs. (33)- (34) in Borrero (2007). These equa- 
tions take into account how the external field bends and 
surrounds the flux tube. Although the model is three- 
dimensional, variations along the x-axis are neglected. 
The model is described by 6 parameters: Bq and 70 
(strength and inclination of the external magnetic field 
far away from the flux tube), R and zq (radius and cen- 
tral position of the flux tube), and finally B^to and Vxto 
(component of the flux tube magnetic field along the flux 
tube's axis and magnitude of the Evershed flow). Figure 1 
(top panel) illustrates the magnetic field lines in the plane 
perpendicular to the tube's axis. 

The velocity and magnetic field vectors enter the static 
momentum equation. This yields the density p(?/,z), and 
gas pressure Pg{y^z)^ that ensure force balance (see Bor- 
rero 2007 for details). Once gas pressure and density are 
known, the temperature T{y^ z) is evaluated through the 
equation of state for ideal gases. A varying molecular 
weight is used to account for the partial ionization of the 
different atomic species. 

Once all the relevant quantities are known in the LRF, 
<S, we project the magnetic field and velocity vectors on 
the observer's reference frame (ORE), S" . This is accom- 
plished by a rotation of angle ^ along the vertical z-axis. 
This rotation ensures that the resulting x -axis meets the 
line of symmetry of the sunspot. A second rotation of an- 
gle O along Cy is finally needed to direct the resulting z 
along the observer's line-of-sight. Mathematically, 

5' = {ex^ey^e;'}=7^^.(e) x7^,(^) x5, (1) 
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Fig. 1. — Top panel: magnetic field lines tranversal to the flux 
tube's axis. Note that inside the flux tube the magnetic fleld is 
mostly aligned its axis and therefore it is not shown here. Bottom 
panel: inclination of the magnetic fleld in the observer's reference 
frame (i.e.: inclination respect to the line-of-sight). The white arrow 
represents a possible line-of-sight of an observer at B = 45° looking 
at a penumbral flux tube located at ^ = 45° with respect to the 
line of symmetry in the center side. In this example, the intersec- 
tion of the ray-path with the uppermost point in the yz plane occurs 
at (yo, ^macc)=( — 150, 200) [km]. This point is indicated by a black 
asterisk. The angle between the vertical z-axis (white dashed line) 
and the projected line-of-sight is [3 ^ 35°. Other ray-paths would 
be parallel to the one drawn but intersect at different yo's. 



where lZy'{Q) and lZz{^) are the corresponding rotation 
matrices. In this manner, we can locate the flux tube 
at any azimuthal position ^ within the sunspot, and to 
place the sunspot at any heliocentric angle, O, on the so- 
lar disk. Thus, we can determine the line-of-sight velocity 
'^los = -> magnetic field strength B (calculated in any 
reference frame) , inclination of the magnetic field with re- 
spect to the line-of-sight 7 = cos~^[5^/5], and the angle 
of the magnetic field vector in the plane perpendicular to 
the line-of-sight ^ = tan~^ [5^ l^y\' example of 7(7/, z) 
is presented in Fig. 1 (bottom panel). 

Finally, our physical parameters are expressed as func- 
tion of (2/, 2;), but for our radiative transfer calculations 
we need to know them along the line-of-sight. To this end, 
we project the observer's line-of-sight onto the yz plane: 
lA = sin ^ sin Ocy — cos Ocz. With this vector we can con- 
struct the parametric equation of the line-of-sight in the 
yz plane: 



{y,z) = {yo.Zmax) + A (sin ^ sine, -COS e) (2) 

where (7/0, Zmax) is the uppermost intersection point of the 
ray-path with the yz plane. The slope-intercept form of 
the line-of-sight is: 



y = yo-^ {zmax - z) tan e sin ^ 



(3) 



Note that, either at disk center (6 = 0) or along 
to the line of symmetry of the sunspot (^ = 0,7r/2), 
the ray-paths are given hy y = y^. At any other posi- 
tion, the ray-paths enter the yz plane, forming an angle 
[3 = tan~^(tan sin ^) with the vertical direction. As ex- 
pected, P = ±Q perpendicularly to the line of symmetry 
( ^ = 7r/2,37r/2). An example of the projection of the 
line-of-sight onto the yz plane is presented in Fig. 1. 

After these geometrical considerations, we are now 
ready to solve the radiative transfer equation and ob- 
tain theoretical Stokes profiles from our embedded fiux 
tube model. We have employed three different numerical 
codes: SIR (Ruiz Cobo & Del Toro Iniesta 1992), DIA- 
MAG (Grossmann-Doerth 1994) and SPINOR (Frutiger 
2000) and verified that the results are consistent among 
them. We first calculate the emerging polarization pro- 
files for 128 rays that enter into the yz plane with diff'erent 
?/o's (see Fig. 1). We then compute the Net Circular Po- 
larization (NCP) as the wavelength integral of Stokes V. 
The total NCP is obtained as sum of the NCP produced 
by each individual ray-path (see example in Fig. 2; top 
panel) . 

3. OBSERVATIONS 

Five sunspots at five different positions on the solar 
disk were observed using the Tenerife Infrared Polarimeter 
(TIP; Martinez Pillet et al. 1999) at the Vacuum Tower 
Telescope (Izaha Observatory, Spain), and the Advanced 
Stokes Polarimeter (ASP; Elmore et al. 1992) at the Dunn 
Solar Telescope (Sacramento Peak, USA). The full Stokes 
vector of Fe I 15648.5 A (TIP) and Fe I 6302.5 A (ASP) 
were recorded . For each sunspot, radially averaged az- 
imuthal variations of the NCP, A/'(^), are computed. Here 
^ runs counter-clockwise, with ^ = referring to the line 
of symmetry on the center side of the penumbra. The de- 
termination of the position of the line of symmetry from 
the observations is a difficult task and can only be done to 
an accuracy of about ±10°. 



Table 1 

Summary of observed sunspots. References: a 
(Cabrera Solana et al. 2006), b (Borrero et al. 
2006), c (Bellot Rubio et al. 2004), d (Borrero et 
AL. 2005). 
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Fig. 2. — Top panel: Net Circular Polarization produced by each 
ray-path piercing the flux tube at different y positions. The flux 
tube is located at the line of symmetry on the center side (dashed, 
^ = 0) and on the limb side (solid; ^ = tt). Bottom panel: variation 
of the total NCP (integral over y in the top panel) as a function of 
the heliocentric position of the sunspot: center-to-limb- variation. 

4. OBSERVED VS THEORETICAL NCP 

4.1. Center-to-limb variation 

We have computed the NCP that emerges from our flux 
tube model for different hehocentric angles 0. This is the 
so-called center-to-limb variation of the NCP. We have re- 
peated this calculation for flux tubes located at the line of 
symmetry of the sunspot, on the center side = 0) and 
on the limb side = tt). The NCP was evaluated for the 
same spectral lines of our observations (Sect. 3). Results 
are presented in Fig. 2 (bottom panel). The parameters 
used for the flux tube model are: Bq = B^m = 1000 Gauss, 
7o = 60°, Vxm = 6 km s~^, R = 75 Km, = km. These 
values are consistent with result from spectropolarimetric 
observations of the penumbral flne structure. In addition, 
we have used the hot umbral model by Collados et al. 
(1994) to represent the thermodynamic parameters of the 
atmosphere surrounding the flux tube. 

Martinez Pillet (2000) has presented the observed 
center-to-limb variation of the NCP in Fe I 6302.5 A for a 
large number of sunspots. Our predictions (Fig. 2; bottom 
panel; blue lines) agree very well with his flndings. Our 
theoretical JV{Q) curve for 6302.5 A does not cross zero at 
cos 9 = [0.8, 1], however. This is due to our model simpli- 
flcations: flux tube is always perpendicular to the vertical 
z— axis, and the external atmosphere does not harbor any 
flows. 

For the near infrared neutral iron line at 15648.5 A we 
are not aware of any similar work to that of Martinez Pil- 



let. However, our theoretical curve is in very good agree- 
ment with other theoretical predictions that use a simpler 
uncombed scenario (see Miiller et al. 2002, 2006). The 
reason is that, along the line of symmetry (^ = 0, tt), ray- 
paths are not inclined on the yz plane regardless of the 
heliocentric angle (Eq. 3). 

4.2. Azimuthal variations 

In this section we compare the observed A/'(^), for dif- 
ferent sunspots at different heliocentric angles O, with 
the theoretical predictions from the embedded flux tube 
model. This is done individually for each sunspot. The- 
oretical curves have been obtained using the same model 
parameters as in Sect. 4.1. Results are presented in Figure 
3. It can be seen that our model is able to reproduce many 
features of the observed azimuthal variations of the NCP. 
This achievement is especially remarkable if we consider 
that all we have done is changing the spectral line and 
heliocentric angle (model parameters were kept constant). 

Particularly interesting is the fact that we can also pre- 
dict the existence of secondary maxima/minima in the 
A/'(^) curves. This is clearly the case of Fe I 15648.5 A , 
where more simple uncombed models, predict only the ex- 
istence of 2 maxima and 2 minima at all heliocentric an- 
gles. The improved agreement between observations and 
predictions is to be ascribed to our more realistic model- 
ing of the external magnetic fleld bending and wrapping 
around the horizontal flux tube. This is a crucial ingredi- 
ent for Fe I 15648.5 A , where gradients in azimuthal angle 
of the magnetic fleld play a major role (Landolfl & Landi 
degl'Innocenti 1996; Miiller et al. 2002). 

An important detail is that, in our model, the Ev- 
er shed flow is channeled along the horizontal flux tube. 
This means that the line of sight velocity, vios = = 
— Va^^o cos ^ sin 6, vanishes always perpendicularly to the 
line of symmetry (^ = 7r/2,37r/2). Thus, gradients in vios 
do not exist there, yielding always zero Net Circular Po- 
larization. In agreement with Miiller (2001), decreasing 
the magnitude of the Evershed flow decreases the amount 
of NCP roughly linearly. In the absence of a complete pa- 
rameter study, this effect alone cannot be used to rule out 
weaker horizontal flows, however. 

Borrero (2007) expressed concerns about flux tubes with 
circular cross sections having very smooth variations in 7. 
He pointed out that those variations might be too small 
to generate enough NCP (Sanchez Almeida & Lites 1992). 
Results presented in this work prove those concerns to be 
unfounded. 
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Fig. 3. — Comparison between the observed (solid circles) ,and 
predicted (red solid lines) azimuthal variations of the Net Circular 
Polarization, for two sunspots at heliocentric angles: B = 27°, 38°. 
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Fig. 4. — Same as Fig. 3 but for three sunspots at = 
40°, 50°, 60°. 

5. CONCLUSIONS 

We have employed the model of Borrero (2007) to pre- 
dict the behavior of the Net Circular Polarization in the 
sunspot penumbra. This model finds the equilibrium con- 
figuration of a horizontal fiux tube with circular cross sec- 
tion, that carries the Evershed fiow, and is embedded in 
a atmosphere with a potential magnetic field pointing to- 
wards a different direction. Energy transfer is neglected 
and therefore we cannot address how the penumbra is 
heated. This model consistently considers how the exter- 
nal magnetic field opens and bends in order to accommo- 
date the horizontal fiux tube. This generalizes the work of 
Solanki & Montavon (1993), Martinez Pillet (2000), Miiller 
et al. (2002; 2006) and Schlichenmaier et al. (2002), by 
considering the 3D geometry of the problem. 

We have compared our predictions with the observed 
NCP in two neutral iron lines and in five different sunspots. 
The agreement between theory and observations is remark- 
able, improving previous determinations based on simpler 
realizations of the uncombed penumbral model. Other 
models for the penumbral fine structure (Sanchez Almeida 
2005; Scharmer & Spruit 2006) should also try to explain 
these observations. In the future, we will attempt to re- 
produce also the full polarization profiles of these spectral 
lines (cf. Borrero et al. 2005,2006; Bellot Rubio et al. 
2004). 

We wish to thank Daniel Cabrera Solana for kindly pro- 
viding some of the observations used in this work. 
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